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ABSTRACT 

CLASS (Conmunl cations Link Analysis and Simulation System) Is a software package developed for 
NASA to predict the communication and tracking performance of the Tracking and Data Relay Satellite 
System (TDRSS) services. This paper describes some of the modeling techniques used in CLASS. 

1. INTRODUCTION 


The Communications Link Analysis and Simulation System (CLASS) presently under development for 
NASA/Goddard Space Flight Center (GSFC) is an Integrated set of FORTRAN programs capable of predicting 
the compatibility and performance of the communication and tracking links for all services and signal 
formats supported by the Tracking and Data Relay Satellite System (TDRSS) . CLASS contains detailed 
models for the TDRSS spacecraft and ground terminal hardware. Including Shuttle-unique equipment, and 
models for the effects of the transmission medium (rain attenuation, multipath, radio frequency 
Interference (RFI). CLASS allows the modeling of the transponder of a TDRSS user such as Shuttle 
either based on a specific hardware Implementation or based on a set of parameters describing the 
signal characteristics at the RF interface between TDRSS and the user transponder. These parameters 
Are used by NASA/ GSFC to specify the quality of the user’s signal at this Interface and are, therefore, 
referred to as user constraints [1]. 

CLASS is capable of verifying the compatibility of a particular transponder design with the TDRSS 
signal formats and predicting the system performance in terms of all the performance parameters of 
Interest to the user. These include data integrity (bit error probability, bit slippage probability, 
probability of carrier phase slips In MPSK systems), synchronization (tracking jitter, slip rates, loss 
of lock probability), tracking (range and range rate), and acquisition (acquisition time, probability 
of false lock) performance. 

The capabilities and structure of CLASS are presented In more detail In another paper In these 
proceedings , Communications Link Analysis and Simulation System" by Robert Godfrey. The effort made 
to validate CLASS Is described In another paper In the Proceedings, "Validation of the Communications 
Link Analysis and Simulation System (CLASS)" by the same authors as the current paper. 

The purpose of this paper Is to describe some of the modeling techniques used In CLASS. The 
components of TDRSS and the performance parameters to be computed by CLASS are too diverse to permit 
the use of a single technique to evaluate all performance measures. Hence, each CLASS module applies 
the modeling approach best suited for a particular subsystem and/or performance parameter in terms of 
model accuracy and computational speed. It was one of the challenges of the CLASS development to 
design a software structure which allows these diverse modules to share one system database. 

■ 

The following sections provide a brief description of the modeling techniques used for four major 
parts of CLASS: the bit error rate performance computation (Section 2), the synchronization/tracking 

subsystem (Section 3), the acquisition subsystem (Section 4), and the evaluation of RFI effects 
(Section 5). 


2. COMPUTATION OF BIT ERROR RATE 


The channel model for the bit error rate analysis must account for the signal distortion occurring 
In the transmitter, the relay satellite, as well as the receiver. A typical return link channel model 
Is shown In Figure 1. The complexity of this channel model, particularly the mixing of linear and 
nonlinear elements, prevents the use of strictly analytical performance c evaluation techniques 
Similarly, the low bit error probabilities of Interest (the design BER Is 10~ 5 > make Monte Carlo type 
simulations prohibitively slow. The approach used In CLASS combines elements of both techniques. The 
signal Is represented as a sampled waveform which allows the modeling of all user constraints? Linear 
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signal distortions are easily incorporated through a fast Fourier transformation (FFT) of the waveform 
and appropriate processing in the frequency domain. The interaction between the signal, adjacent chan- 
nel waveforms, and thermal noise in nonlinear elements (power amplifiers, limiters) Is modeled through 
a modified Barrett-Lampard expansion developed for this program and described elsewhere [2 J . 

The receiver model operates on' the sampled signal and a statistical representation of the thermal 
noise and interference. First, the effect of the signal characteristics on the synchronization systems 
(PN code, carrier, clock recovery) is computed. The effect of these subsystems on the recovered 
symbols is then Included in the bit error performance evaluation by computing first the BER conditioned 
on the various synchronization errors and then averaging over the appropriate probability densities. 
The BER in the case of coding and 3-level quantization at decoder Input is computed via the Rg 
approximation [3]. The channel cut-off rate Rg, computed from the channel output probabilities 
conditioned on channel Input, is assumed to fully characterize the channel. The BER is assumed to be a 
function of code rate and Rg only, those functions being given In [*]. 

The above approach Is not intended as a universal modeling tool for nonlinear channels. The 
modified Barrett-Lampard expansion describes the spectral characteristics of the output of the 
nonlinearity and it requires an excessive number of terms when highly nonlinear elements, such as hard 
limiters, are Included in the channel model. Similarly, the approximations made in the statistical 
representation of the noise and Interference become inaccurate when highly nonlinear elements are 
modeled. However, for the characteristics of the TDRSS channels the accuracy of the models was 
verified by analysis and comparison to Monte Carlo simulations. 

To point out the Importance of a comprehensive link model, the single parameter sensitivities with 
a perfect signal and a linear, wideband channel are'compared In Figure 2 with sensitivity results 
obtained with a typical TDRSS link model. .The plots show that single parameter sensitivity results 
with a perfect signal and a linear, wideband channel can be quite misleading. Similarly, Table 1 shows 
that for another typical link the sum of the degradations that result from increasing each of several 
parameter values separately can be different from the degradation that results from increasing all the 
values at once. 


3. MODELING OF THE SYNCHRONIZATION/TRACKING OPERATION 

The synchronization subsystems modeled are the PN code, carrier, and clock recovery loops In their 
various Implementations In the TDRSS ground station. Several performance parameters are of Interest: 

(1) the tracking jitter due to thermal noise and signal amplitude and phase variations, 

(2) the tracking offset due to loop stress and vehicle dynamics (Doppler, acceleration), 

(3) the rate of cycle slips due to all the above sources and to untracked clock jitter, 

(4) the probability of dropping lock due to vehicle and antenna dynamics and signal dropouts. 

These subsystems are modeled by a combination of analysis and simulation. The analysis accounts 
for the effect of thermal noise, phase noise (clock jitter), and vehicle and antenna dynamics on the 
synchronization performance (based on nonlinear tracking models [5]), while the signal waveform 
generated for the BER performance evaluation Is used to characterize the effect of waveform 
distortions. The phase- tracking jitter generated by these signal distortions Is generally referred to 
as pattern Jitter. By linearly combining It with the thermal noise effect the overall subsystem 
performance Is obtained. 

The characterization of the tracking performance requires. In addition to the above aspects, an 
accurate model of the low-frequency phase noise In the TDRSS and user transponder. All TDRSS links can 
operate in a coherent turnaround mode. This means that the phase noise processes on the ground. In the 
TORS, and in the user equipment are correlated. However, due to the small phase error values of 
Interest, linear models apply and the system can be characterized by a linear network. A thorough 
description Is provided In [6J. 


4. MODELING OF SIGNAL ACQUISITION 


The signal acquisition process includes the search for the PN code epoch (for spread spectrum 
signals), the carrier frequency and phase acquisition, and the clock synchronization. In addition, the 
Ku-band links require the spatial acquisition of the autotrack system. 

Analytical models were developed for the different hardware Implementations used in the TDRSS 
ground station and various user transponders. The signal model used reflects the user constraints as 
well as the linear and nonlinear distortion effects of the channel. 

The PN acquisition model computes the characteristic function of the search time based on the a 
priori code epoch uncertainty and the search algorithm. Specific cases of interest Include single 
dwell time, dual dwell time, and variable dwell time systems for a circular search over the whole code 



or an expanding window search over a small part of the long code (e.g. 
detailed description of the analytical approach Is given In [7,8]. 


during reacquisition) . A 


The carrier acquisition algorithms modeled Include frequency-locked loop (FLL) aiding, swept 
acquisition, and frequency pull-in (self acquisition). The modeling Is based on analytical models 
[5,9,10] with appropriate modifications to account for signal distortions (user constraints and channel 
effects). The FLL may operate with a low SNR at the Input to the frequency discriminator. The model 
includes, therefore, a spike noise component at the discriminator output. The FLL bandwidth Is small 
however, which permits the use of a linearized loop model. 


A major component of these models is the acquisition/tracking monitor and the transition between 
the two modes. The response time of the monitor can contribute significantly to the overall 
acquisition time. Also, the Implementation of this monitor determines the probability of false lock. 


The clock recovery loops modeled are all self-synchronizing, 
these systems do not contribute significantly to the acquisition time. 


With typical clock stabilities. 


5. MODELING OF RFI EFFECTS 

It Is expected that some of the TDRSS links will be subject to high-powered RFI. This 

interference may be a mixture of pulsed CW signals and noise pulses of a few microseconds duration with 
random arrival times. t 

CLASS Is designed to account for the RFI effect on all aspects of TDRSS performance. First we 
discuss modeling of the effect on BER, then on synchronization and acquisition, and last on tracking. 

Two approaches were taken In the case of BER. First an analytical model was developed based on 
assumptions that hold for the general TDRSS S-band user. The model breaks down for the Shuttle S-band 
return link In the prediction of decoder performance. For this reason, a second approach, a Monte 

Carlo type simulation, was also taken. A benefit of having both approaches Is that they serve to 
verify each other. 


JH e ^ na1yt1ca1 a PP ro , ach 1n the case of BER Is based on the sample-sum approximation to the 
matched-filter output [11]. The approach is implemented In different ways depending on whether the 
symbol rate Is high or low relative to the Inverse of the RFI pulse duration. In both cases the 
characteristic function of the matched-filter output is obtained. From It, In the case of a 
convolutlonally encoded link, the cut-off rate Rq Is computed, from which In turn BER Is computed. An 
early account of the modeling approach Is given In [12]. However, extensive testing demonstrated the 
need for a host of refinements to Insure good accuracy for all RFI environments and data rates. These 
updated models are described below. 


In the high-rate analytical model for BER at most one RFI pulse Is assumed to occur In a symbol 
duration. The Gaussian or non-Gausslan characteristic function of the matched-filter output Is 
obtained, conditioned on no RFI or on the type and power of RFI occurring In the symbol. For CW RFI 
conditioning Is also on the RFI frequency and Initial phase. The conditional characteristic functions 
are weighted and summed to yield the unconditional function. Since Interleavers are used by the 
general TDRSS S-band user (but not Shuttle) for high symbol rates. It can be assumed for the decoder 
model in that case that the channel is memoryless. 

In the low-rate analytical model for BER, rich RFI pulse Is assumed to be wholly contained In one 
symbol. The characteristic function of the contribution to the matched-filter output of one RFI pulse 
Is computed, conditioned as above. This characteristic function Is normalized by the characteristic 
function of no RFI occurring for the pulse duration. The average of all such functions Is obtained. 
Poisson pulse arrival statistics are employed to form the non-Gausslan characteristic function of the 
matched-filter output as the product of a term Involving the function just described and the no-RFI 
characteristic function. 


In the Monte Carlo type simulation to model the RFI effect on BER, no Interleaver or deinterleaver 
Is present, which reflects the case of the Shuttle return links. 

The models of the RFI effect on all the synchronization subsystems and most of the acquisition 
subsystems are similar. In these subsystems there are a filter with bandwidth roughly equal to the 
symbol rate, a nonlinear element and possibly some other elements, and then a filter with a bandwidth 
less than the lowest symbol rate. As In the BER model, two analytical models are actually used. For 
low symbol rates the first filter Is assumed to average the RFI conditions. For high symbol rates the 
first filter Is assumed wide enough to pass the RFI undistorted and the second filter Is assumed narrow 
enough to do the averaging. The results are Insensitive to the exact symbol rate that Is chosen for 
the cross-over point between the models. For swept carrier acquisition the time during which the loop 


_800 



can pull In Is assumed long enough to contain many pulses of each type of RF1. 

The RFI affects Doppler tracking by Increasing carrier phase Jitter and cycle slips. 
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Figure 1. TYPICAL LINK MODEL. 
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Figure 2a 


Figure 2b. 


Figure 2c. 
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COMPARISON OF PARAMETER SENSITIVITIES FOR CLASS CHANNEL MODEL AND IDEALIZED CHANNEL 
MODEL FOR MODULATOR PHASE ERROR. 



COMPARISON OF PARAMETER SENSITIVITIES FOR CLASS CHANNEL MODEL AND IDEALIZED CHANNEL 
MODEL FOR RECOVERED CLOCK PHASE ERROR. 
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Table 1. CNR LOSS DUE TO SEVERAL PARAMETERS EXCEEDING NOMINAL VALUES (COURTESY OF 
ROBERT GODFREY OF NASA). 




Parameter 

Nominal Value 

Actual Value 

CNR Loss 


1 . 

Modulator Phase Imbalance 

3° 

7° 

0.2 dB 

■ '• .v'V' 1 

— - v:. *' i 

2. 

Modulator Gain Imbalance 

+0.25 dB 

+0.8 dB 

0.87 


3. 

Phase Nonlinearity 

+3° 

+9° 

0.5 


4. 

Gain Flatness 

+0.3 dB 

+3 dB 

— 


5. 

Gain Slope 

+0.1 dB/MHz 

1 

0.15 


6. 

AM/PM 

1 2°/dB 

20°/dB 

0.58 


7. 

Phase Noise 

1° 

3° 

0.1 



Composite Impact 

i 


3.18 dB 


(BPSK, NRZ, Rate-1 /2 Coding, No PN Spreading) 


> 
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